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ABSTRACT 



We report ionization cross section measurements for electron impact single 
ionization (EISI) of Fe^^"*" forming Fe^^"*" and electron impact double ionization 
(EIDI) of Fe^^"*" forming Fe^^"*". The measurements cover the center-of-mass en- 
ergy range from approximately 230 eV to 2300 eV. The experiment was per- 
formed using the heavy ion storage ring TSR located at the Max-Planck-Institut 
fiir Kernphysik in Heidelberg, Germany. The storage ring approach allows nearly 
all metastable levels to relax to the ground state before data collection begins. 
We find that the cross section for single ionization is 30% smaller than was previ- 
ously measured in a single pass experiment using an ion beam with an unknown 
metastable fraction. We also find some significant differences between our ex- 
perimental cross section for single ionization and recent distorted wave (DW) 
calculations. The DW Maxwellian EISI rate coefficient for Fe^^+ forming Fe^^"*" 
may be underestimated by as much as 25% at temperatures for which Fe^^+ is 
abundant in coUisional ionization equilibrium. This is likely due to the absence 
of 3s excitation-autoionization (EA) in the calculations. However, a precise mea- 
surement of the cross section due to this EA channel was not possible because 
this process is not distinguishable experimentally from electron impact excita- 
tion of an n = 3 electron to levels of n > 44 followed by field ionization in the 
charge state analyzer after the interaction region. Our experimental results also 
indicate that the double ionization cross section is dominated by the indirect 
process in which direct single ionization of an inner shell 21 electron is followed 
by autoionization resulting in a net double ionization. 



- 3 - 



Introduction 



CoUisionally ionized atomic plasmas are formed in a range of astrophysical objects 
including stellar coronae, galaxies, and supernova remnants. The charge state distribution 
(CSD) of such plasmas is determined by a balance between electron impact ionization 
(EII) and electron-ion recombination. The CSD plays an important role in a wide range 



of spectroscopic diagn ostics used 



:o infer electron temperature , electron density, a nd 



elemental abundances dBrickhouse 



1996 



Landi fc Landini 



1999 



Bryans et al. 



20091 ). Most 



CSD calculations are carried out under the assumption of collisional ionization equilibrium 
(CIE) for conditions where the electron density is low (meaning three body recombination 
is unimportant), radiation can be ignored, and there is no dust. In this case the ionization 
and recombination rates are the same so that neUqaj = nen^+ia^^^, where rie is the electron 
density, Ug is the density of ions of a particular element with charge q, al is the rate 
coefficient for ionization from g to g + 1, and is the rate coefficient for recombination 
from g + 1 to g. Rewriting gives nq^i/nq = al/a^^^, thus making clear the importance of 
accurate ionization and recombination data to model the CSD of a plasma. 

For CIE, generally only electron impact single ionization (EISI) need be considered since 
for a given charge state multiple ionization is signi ficant only at temperatures so high that 



the fractional abundance of that ion is negligible (ITendler et al. 



19841 1. Multiple-electron 



ionization, such as electron impact double ionization (EIDI), need be considered only when 
modeling the CSD in dyna mic systems where ions are suddenly exposed to higher electron 



temperatures (IMii. 



(IReale fc Orlando 



lei 



19861 ). Examples of such non-equilibriura 



2008|). supernova remnants (jPatnaude et al. 



syste ms include solar flares 



clusters flAkahori fc Yoshikawa 



2010). 



20091 ) ■ and merging galaxy 



by 



Much wor k has been carried out in deriving the necessary EII data (see the compilation 



Dere 



20071 ). however sizable discrepencies in the data exist even among the most 
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recent compilations (iBryans et al.l 120091 ) . Experimental EII measurements, like those 
reported here, will help to resolve these differences. These measurements can also be 
used to benchmark theory enabling more accurate EII cross section calculations for ions 
isoelectronic to those measured. 

The difficulty of producing well characterized ground state ion beams has been 
a major limitation for obtaining the accurate EII data needed for CIE models. Most 
measurements of EII have been performed in a single-pass geometry in which metastable 
ions generally have not had enough time to radiatively relax to the ground state before the 
EII measurements are performed. In this paper we describe an experiment employing an 
ion storage ring. This arrangement allows the ions to be stored long enough for typically 



all the metastable levels to r adiatively relax be 



been exploited previously by 



Linkemann et al 



bre da ta acquisition begins. This fact has 



fll995a 



Here we study ionization of Fe This ion produces strong spectral lines which can 



be used for plasma diagnostics and instrument calibration (jPel Zanna fc Mason 



is also a part 



Brown et al. 



icularly important ion for observations of the solar corona (IMoses et al 



2QQH\. It 



1997 



20081). 



e: 



SI of Fe has already been the subject of some theoretica 



work (I Younger 



1983 



Pindzola et al. 



1986 



Gregory et al 



1987 



Dere 



and e xperimental 



20071). 



Gregory et al 



(Il987l ) measured the ionization cross section for Fe^^"'' using a single pass experimental 
geometry. The EISI cross section measured by that experiment was about 30 % larger tha n 



that predictec 



Gregory et al. 



)y th e distorted wave calculations of 



Pindzola et al 



fll986h or 



Derel fl2007h . 



(119871 ) compared their experimental result to theoretical calculations for 
ionization from the ground and metastable levels and concluded that the discrepancy was 
due to a large metastable population in the ion beam in the experiment. The measurements 



reported here, by being able to generate nearly pure ground state ion beams, can resolve 
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this issue and also provide an excellent test case to compare crossed beams experimental 
results with the storage ring technique. 

We study in detail EISI of P-like Fe^^^ forming Si-like Fe^^^ over the electron-ion 
collision energy range of 230 — 2300 eV. This energy range includes the following direct 
ionization channels: 

Fei2+(2s2 2/3s3p3) + 2e- 

. (1 

Fei2+(2s2 2/ 3s2 3^^) + 2e- 
Fei2+(2si 2/ 3^3^ ^ 2e- 



e" + Fe^^+(2s^ 2^** 3s^ 3/) ^ <( 



The energy threshold s for direct ionization a re 330.79 eV from the 3p subshell and 357.40 eV 



from the 3s subshell (IRalchenko et al. 



I2OO8I ). Excitation-autoionization (EA) can also occur 
starting at the ionization threshold of 330.79 eV through the electron impact excitation 
of a 3s electron to an autoionizing state. We are unaware of any theoretical work on this 



specific EA channel for Fe^^"*". 

Theoretical calculations have shown that EA is expected to make an important 
contribution to the ionization cross section above the excitation threshold for n = 2 



electrons at ~ 710 eV (IPindzola et al. 



19861 ). In that same energy range resonant processes 
such as resonant excitation double autoionization (REDA) and resonant excitation 
auto-double-ionization (READ!) are predicted in the cross section. These occur through 



dielectronic ca pture when an ion forms an excited state that 



two electron s (ILaGattuta fc Hahn 



1981 



Miiller 



Henry fc Msezane 



subsequently decays by e 



1982 



Miiller et al 



1988a b 



ecting 



20081 ). In the REDA process the electrons are released sequentially, whereas in 



READI the two electrons are ejected simultaneously. 

The theoretical energy thresholds for ionizatio n of the 2p or 2s electrons are 1073 eV 



and 1199 eV, respectively (IKaastra fc Mewe 



19931 ). Ionization from the 2p subshell 
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radiative 



y stabilizes with a pred icted probability of 1.1% and from the 2s with a probability 



19931 ). In both cases it is more likely that direct ionization of 



of 8.0% flKaastra fc Mewd 
an n = 2 electron will be followed by autoionization giving a net double ionization to form 
Fe^^"*". Thus, ionization of an inner shell electron is expected to be a very small contribution 



to the total sing . 



f Miiller &: Frodl 



e ion ization cross section, but it can be significant for double ionization 



19801). 



Electron impact double ionization of Fe^^"*" forming Al-like Fe^'^"^ was measured over the 
same energy r ange of 230 - 2300 eV. The energy threshold for EIDI through direct ionization 



is 691.83 eV (IRalchenko et al. 



20081 ). As discussed above, double ionization can also occur 
through single ionization of an inner shell electron with subsequent autoionization. Other 
indirect processes such as EA and resonant processes analogous to those expected fo r single 



ionization may also occur in the double ionization cross section (iMiiller et al 



1988a). 



The rest of this paper is organized as follows: In Section [2] we describe the experimental 
setup. The data analysis is presented in Section [3] and uncertainties in Section HI 
Experimental results and comparison to theory are presented for single ionization in 
Section [5] and for double ionization in Section [61 A summary is given in Section [71 



2. Experimental Setup 



EII measurements were performed at the TSR heavy ion storage ring of the 



Max-Plank-Institut fiir Kernphysik in Heidelberg 



, Germany. 



followed the procedure of p r evious experiments (jKilgus et al 



1995a 



Kenntner et al 



1995 



Schippers et al. 



2001 



he experiments basica. 



1992 



Lestinsky et al 



specifically to ionization are described in iHahn et al 
aspects relevant to the present work. 



Linkemann et al. 



ly 



20091). Details related 



(120 10[ ). Here we describe additional 
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The TSR facility is equipped with two separate electron-ion merged beams sections. 



For the resu 


ts presented here the 


(Steck et al. 


1990; 


Pastuszka 


1996) 



studies at tunable rel ative energies. The second electron beam device is known as the Target 
2004! ). This device was oper ated as an e lectron cooler with the energy fixed 



( ISprenger et al 



at what is referred to as the cooling energy ( iPoth 



19901 ). Fixing the Target electron beam 



energy allowed the ion beam to be cooled continuously. This inhibits expansion of the beam 
due to warming and counteracts the drag force, both of which arise as the electron energy 
of the probe beam is varied during measurement. 

A beam of 150 MeV 56pgii+ jQ^g ^g^g injected into TSR in a series of 5 pulses spaced 
0.8 s apart. The ion beam was merged with the two electron beams described above. 
Initially both the Cooler and Target electron beams were set to the space-charge-corrected 
cooling energy of 1460 eV. The center-of-mass energy spread is limited by the Cooler 
electron temperature, which can be described by a flattened Maxwellian distribution with 
temperatures in the perpendicular and parallel directions of typically k^T^ = 13.5 meV and 
k^Tw = 180 fieV (Kilgus et al. 1992; Schippers et al. 2001; Novotny et al., in preparation). 
During measurement, the average ion current was 1-2 fiA. 

A delay of 2 - 3 s followed the last injection pulse before data collection began. This 
allowed Fe^^"*" metastable states to decay. To estimate the metastable fraction during 
measurement we modelled the level populations as a function of storage time starting from 
a Boltzmann distribution with temperature 750 eV, corresponding to the approximate 
collision energy of the electrons in the stripping foil as th e ions pass through the foil in the 



accelerator. The energy levels and Einstein coefficients of 



Ralchenko et al. 



(120081 ) were used 



for the calculation. The longest lived metastable le vel was the ^D^/ 



configuration, which has a lifetime of about 0.5 s ( iRalchenko et al. 



leve 



within the ground 



20081 ). The metastable 
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fraction of the ion beam at the end of the initial coohng cycle is estimated to be < 0.5%. 

Products of ionizing electron-ion collisions formed in the Cooler were separated from 
the Fe^^"*" parent beam by the first dipole magnet downstream of the interaction section. 
Recombination and ionization were measured simultaneously using separate detectors on 
opposite sides of the parent ion beam. The detectors were positioned for maximum signal 
collection by stepping each unit horizontally and vertically in small increments across 
each product beam. Single or double ionization events were detected by positioning the 
ionization detector to intercept either the Fe^^"*" or Fe^'^"^ product beam. 



The ionization dete ctor employs a channel electron multiplier (CEM; 



Rinn et al. 



1982 



Linkemann et al. 



1995bl ). A suitable CEM discriminator level was determined by measuring 
the total pulse height distribution. The level was set so that the detection efficiency 
of ionization events was essentially unity. The total ionization detector count rate was 
< 20 kHz with a typical pulse width of 30 ns. As a result, the signal loss induced by the 
dead time of the CEM was negligible. 

Data acquisition began after the injection and initial cooling cycle already described. 
During data collection the relative energy between the probe electron beam and the ion 
beam was varied. Each energy scan consists of ~ 250 - 700 pairs of steps, one step at the 
measurement energy and the other at a fixed reference energy. The total duration of each 
step was 12-35 ms. There was a delay of 5 - 20 ms at the beginning of each energy step 
before data were collected in order to allow the power supply to settle at the new voltage. 
Data were collected for the remaining 5-15 ms of each step. 

The laboratory electron energies were always chosen to be higher than the cooling 
energy and fell between 4000 and 8000 eV. Each energy scan covered a range of 80 - 2000 V 
in the laboratory frame. A fast high- voltage amplifier with a dynamic range of ±1000 V was 
used to quickly switch voltages for the energy scans. This high-voltage amplifier was used 



- 9 - 



in combination with a slower power supply to lift the fast amplifier into the voltage range 
desired for the energy scan, during which the slower power supply maintained a constant 
voltage. 

The number of pairs of steps in the scan and the range of the energy scan were chosen 
to balance the desired energy resolution against the lifetime of the ion beam, which was 
about 18 s. The energy scans used fell into three broad categories. Low resolution overview 
scans, covering the entire 2000 V range permitted by the fast amplifier, revealed the general 
shape of the cross section. Medium resolution scans, covering an ^ 330 eV laboratory 
energy range, were used to capture major features of the cross section, such as EA, and to 
fill in the details in the overview scan. Finally, high resolution scans, covering a range of 
~ 80 eV in the laboratory frame, were used to resolve even finer details such as resonances. 
The medium and high resolution scans were normalized to the low resolution scans to 
correct for ion current measurement offsets as described in Section HJ 

In between every two measurement steps a reference step was used to estimate the 
background by measuring the count rate at a fixed reference energy. For single ionization 
measurements the reference energy was set below the single ionization threshold for energy 
scans up to a center-of-mass energy of about 950 eV. In this case the single ionization count 
rate at the reference energy was only due to single electron stripping (SES) off the residual 
gas. In the case of double ionization, energy scans extending to 1400 eV were performed 
with the reference below the double ionization threshold so that the count rate at reference 
was due only to double electron stripping (DES). At energies above 950 (1400) eV the 
limited dynamic range of the fast high voltage amplifier prevented us from setting the 
reference point below the single (double) ionization threshold. For these higher energies the 
reference point was set to an energy where the cross section had already been determined 
from lower energy scans. For all data runs the background count rate due to stripping at 
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the reference energy was corrected to better represent the background at the measurement 
energy using the method discussed in the next section. 

Cycles of ion injection, coohng, and energy scan were repeated for ~ 1 hour to improve 
the statistical accuracy. After many cycles the process was repeated for a new energy 
range. Energy ranges were chosen to maintain at least 50% overlap with other scans for 
high statistical accuracy and to correct for systematic offsets due to fluctuations in the ion 
current measurement calibration. 



3. Data Analysis 



Hahn et al. 



(120101 ). To review 



The data analysis follows the procedure described in 
briefly, the cross section for single ionization ai versus energy is obtained from the measured 
ionization rate coefficients {(JiVrei), which are averaged over the velocity spread Afrei of 
the experiment. Because the center-of-mass energy spread is very small, Av^ci "^^rei 
and the cross sections can be calculated by dividing the averaged rate coefficients by the 
relative velocities. The cross section for single ionization, in terms of measured quantities, 
is therefore (see Appendix E]) 



(^l[Em) = < 7Z M r^.rn^T + l^I^rel) l^r. 



V,,i { [1 - P-MEm)]nfNrL ' ' nf[l- A/3c(^rn)] 

Here Rf^{E^) denotes the total single ionization count rate at the measurement step; 
Ri{E^) denotes the background single ionization count rate at the measurement step. 
which is pr o portio nal to the count rate at the reference step -Ri(-E'r) as described in 



(2) 



Hahn et al. 



(120 lOf ): L = 1.5 m is the length of the interaction region for the probe beam; 
the factor (1 — (3i(3e) is a relativistic correction where /3i and /3e are the ion and electron 
velocities normalized by the speed of light, fj/c and fe/c, respectively; and Ej. are the 



center-of-mass energies at measurement and reference, respectively. The electron density 



-li- 



ne is of the order of 10^ cm ^ and is calculated from the measured electron current and 



the geometry of the probe beam (IKilgus et al 



19921 ). The total number of stored ions per 



unit length A^i is calculated from the measured ion current. The ElSl rate at the reference 
energy (o"ifrei) (-Er) is included in order to remove the EISI contribution to the background 
count rate R\ when the reference energy is above the EISI threshold. The factor n\/rif 
accounts for the different electron densities at reference and measurement. The ElDl cross 
section ctdi can be calculated in the same manner as the EISI cross section by replacing the 
single ionization count rates in equation ([2]) with the analogous double ionization rates and 
the term (cxifrei) (-^r) by the equivalent double ionization term (cxDifrei) (-^r)- 

For all measurements, the electron and ion beams also interact in the merging and 
demerging sections on either side of the straight section in which the electron probe and ion 
beam co-propagate. In these sections the ion and electron beams meet at an angle causing 
the relative velocity to be greater than when the beams are colinear. These toroidal effects 



Lampert et al. 



(Il996l ) to correct an overestimate of 



are accounted for using the method of 
the cross section of ~ 20%. 

The total ionization detector signal at the measurement energy is made up of EISI or 
EIDI plus a stripping background. We estimate the background using the rate measured 
at a reference energy where crj or o"di is either zero or known from previous measurements so 
that the stripping rate can be inferred. For EIDI the background is primarily due to DES, 
analogous to the situation in the EISI measurement where the background is due to SES. 
For the double ionization measurement it is possible for the combination of SES plus EISI 
in sequence to contribute to the background, thereby introducing an energy dependence 
that cannot be accounted for by measuring the count rate at the fixed reference energy. 
However, the expected rate for this multiple collision process was estimated based on the 
measured EISI and SES count rates to be at most about 10~^ times the rate of DES. Thus, 
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multiple collisions are extremely rare and can be ignored in the analysis. 

The background rate R\ is proportional to the SES or DES cross section, the number 
of ions, and the residual gas density of the vacuum. The stripping cross sections depend 
only on Vi, which remains constant throughout the measurement. However, the ion beam 
decays during the measurement, so there is a slight difference in ion number related to the 
time delay between the measurement step and the reference step. Additionally, the residual 
gas density can vary systematically with energy in a way that seems to be related to the 
probe beam electron current. This can also introduce a systematic distortion of the cross 
section as a function of energy. Taken all together we expect that R\{E^) ^ R\{Ej.). 

We corr ect for both of the above systematic errors in our analysis using the method 



described in 



Hahn et al. 



( I2OIOI ). The end result of the correction procedure is to adjust the 
shape of the cross section by about 3%. The correction uses the recombination signal at 
high energies as a proxy for the pressure to correct the reference count rate so that it better 
reflects the true background rate at measurement. This works because at high energies 
the recombination signal is dominated by single electron capture (SEC), which depends on 
residual gas density and ion current in the same way that SES and DES do. The electron 
capture signal can therefore be used to detect relative differences in pressure and ion current 
between the reference and measurement steps. The analysis of the recombination data 
shows that dielectronic recombination (DR) of Fe^^^ contributes less than 10% of the total 
recombination signal above the ionization threshold (Novotny et al., in preparation). This 
measured DR component is subtracted from the total recombination signal in order to use 
only the SEC signal for the correction. 

When the reference energy is greater than the ionization threshold the above procedure 
has to be modified. It is possible to account for the decay of the ion beam using direct 
measurements of the ion current. The ion beam decays exponentially with a characteristic 
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decay time measured to be ~ 18 s. The time between the measurement and reference steps 
was about 20 ms. Therefore, the expected difference between the measured reference rate 
and the actual background rate from stripping at the measurement step due to the decay of 
the ion beam is only 0.2%. The effect of the energy dependent pressure variation was not 
corrected for scans where the reference energy is above the ionization threshold. Based on 
the size of the correction applied to the low energy scans, we estimate that not correcting 
single ionization scans above 950 eV or double ionization scans above 1400 eV introduces a 
distortion of about 3% in those energy ranges. 



4. Uncertainties 

The statistical uncertainty on the cross section is not the same for all energies. For 
single ionization the statistical uncertainty is about 10% at 400 eV and drops to 1% by 
700 eV as the number of counts increases with the increasing ionization cross section. The 
statistical uncertainty remains at that level until 1400 eV, where it grows to ~ 3%, because 
at such high energies the cross section is relatively featureless and we performed only a 
few low resolution energy scans. The double ionization cross section was measured with a 
few low resolution scans, and the cross section is much smaller than for single ionization. 
Consequently, the statistical uncertainties for those results are ~ 5% on average, with larger 
uncertainties below 1000 eV where the cross section and corresponding count rates are very 
small. 

There is a systematic error due to the uncertainty in the stored ion current 
measurement. Here the ion current w as measured non-destructively using a beam profile 



monitor (BPM; 



Hochadel et al. 



1994j ). The absolute calibration of this instrument depends 
on the residual gas pressure and any electronic drifts. The calibration drifts with a timescale 
that appears to be one to several hours. To correct for the drift of the BPM calibration we 
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performed an initial calibration of the BPM. This calibration was seen to be stable over a 
three hour period. During this time period we performed a long range energy scan, which 
covered the energy range of 260 - 950 eV. This scan gives an accurate measurement for 
the shape of the cross section as a function of energy. For all other data runs we adjusted 
the ion current calibration in order to produce agreement with this long energy range 
scan. This normalization is possible because all other quantities were measured much more 
accurately than the ion current, so we can attribute discontinuities between the energy 
scans, which overlap in energy range by 50% and so should give identical results, as being 
caused by shifts in the ion current calibration. Thus, the details of the cross section can be 
filled in without introducing distortions to the shape of the curve from the normalization. 
Finally, we quantified the systematic uncertainty in the magnitude of the cross section 
by repeating the analysis of the low resolution scans using high and low estimates for the 
BPM calibration. We estimate the average la systematic uncertainty from the ion current 
measurement to be about 12%. 

Other sources of systematic error were negligibl e compared t o tha t from the ion 



Hahn et al. 



(120101 ) the ion beam was 



current measurement. In the earlier measurements of 
not cooled during measurement and the expansion of the ion beam introduced additional 
systematic errors. These included possible loss of detection efficiency due to ions far from 
the center of the beam not hitting the detectors and uncertainties in the toroidal corrections 
due to the slight differences in the path length through the merging and demerging sections 
of the interaction region experienced by ions at different transverse locations within the 
expanded beam. In the present experiment, beam profile measurements using the BPM 
showed that constant cooling limited the ion beam width to ~ 1 mm. Since the width was 
small compared to the size of the particle detectors, the detection efficiency was essentially 
unity. Similarly, the systematic uncertainty due to the transverse size of the beam in the 
merging and demerging sections was also negligible. The uncertainty on the electron density 
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is about 1% (iKenntnerlllQQSl ). All uncertainties are summarized in Table [TJ 



5. Results for Single Ionization 

The EISI ionization cross section for Fe^^^ forming Fe^^"*" is shown in Figure [H The 
dotted curves in the figure indicate the la systematic uncert ainty due to the ion c urrent 



calibration. The figure also shows the experi mental results of 



to that data used in the CIE calculations of 



Gregory et al. 



Arnaud fc Raymond 



Atomic Code (FAC) distorted wave calculation of 



Derd (120071 1 . 



(119871), the fit 



fll992h . and the Flexible 



Gregory et al. 



(119871 ) is about 



Figure [T] shows that the crossed beams measurement of 
30% larger than the current measurement. This is well outside the uncertainties of either 
experiment, which supports the hypothesis that the earlier measurement had a significant 



metastable population in the ion bear n. Consequently, CSD calcu . 



data are likely to be inacc urate (e.g.. 



Arnaud fc Raymond 



recent CSD c alcula t ion by iBryans et al 
calculation of 



1992 



ations that relied on those 



Mazzotta et al 



19981). A 



( I2OO9I ) used the EISI cross section from the FAC 



Derd (120071 ). Figure [T] shows that this calculation is closer to the present 



experimental results, but significant differences remain. 

In the low energy range from the ionization threshold at approximately 330 eV up 
to about 690 eV, a comparison of our experimental results and theory indicate that the 
ionization cross section is dominated by direct ionization. However, the experimental cross 
section increases faster near threshold than predicted by the calculations. There are two 
related effects that can account for the faster than expected increase. The cross section 
may be enhanced due to electron impact excitation of 3s electrons to states that relax 



by a utoionization. 



ions ([Gregory et al. 



his ionization mechanism has also beeii sugg ested for other P-like 



1983 



Mueller et al 



1985 



Yamada et al 



19881 ). Based on an estimate 
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using the LANL Atomic Code suite (iMagee et al.lll995l ). this mechanism is energetically 
possible for excitation to n > 8 and is sufficient to account for the increased cross section if 
the branching ratio for autoionization compared to radiative stabilization is large enough. 
However, we are unaware of any published data for either EA from 3s excitation or the 
relevant branching ratio. There is also a systematic effect from field ionization that could 
partially account for the discrepancy. The magnetic fields in TSR cause electric fields in the 
rest frame of the ions. If a collision excites an electron to a high enough n level it can be 
field-ionized by the motio nal electric field gener ated in the first dipole magnet downstream 



of the interaction region (jSchippers et al 



200ll ). Here the semi-classical onset for field 



ionization is n > 44. Excitation of a 3p electron to n > 44 can lead to field ionization that 
systematically increases the measured cross section. Excitation of a 3s electron to n > 44 
also leads to field ionization in the experiment, but we expect these states primarily to 
ionize through EA even in the absence of external fields. While the effect of these processes 
on the cross section is most obvious near the ionization threshold, they also increase the 
measured ionization cross section for all higher energies. 

The relative importance of excitation-autoionization versus field ionization can be 
estimated from the dependence of the excitation cross section on the n to which the electron 
is excited. In the Bethe approximation the electron impact excitation cross sectio n for 



dipole allowed transitions from level i to level j is given by (I Van Regemorter 

Stt Jh 



1962|) 



-fijgnao, 



(3) 



where E is the electron collision energy and AEij is the threshold energy for the transition 
both measured in Rydbergs, Jh is one Rydberg, fij is the oscillator strength, g is a. Gaunt 
factor, and is the Bohr radius. For large n the excitation energy AEij is approximately 



constant and equal to the ionization energy. Thus, 



For large n the oscillator strength falls off as 1 / (IBethe fc Salpeted 119571 ) . Therefore the 



or a given colhsion ene rgy aexc,ij oc fij. 
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ratio of the cross section for excitation to field-ionizing levels relative to autoionizing levels 
is approximately (2 Yl'^=44 ^/n^) / {J2'^=s l/""-^) = 0.06. The factor of two roughly accounts 
for the fact that excitation of either a 3s or a 3p electron can lead to field ionization but 
only excitation of 3s can lead to autoionization. This estimate suggests that the systematic 
contribution to the measured ionization cross section from field ionization is small relative 
to excitation autoionization, on the order of a few percent only, but detailed calculations 
taking into account the branching ratios are needed to confirm this. 

From the threshold behavior of the cross section near 700 eV (Figure [2]) we infer that 
EA of n = 2 e l ectron s makes up about 25% of the total cross section at high energies. 



Pindzola et al. 



(119861 ) predict EA from 2p — )■ 3p excitations at 709.3 eV and from 2p — )■ 3d 
excitations at 763.8 eV. This roughly corresponds with the experimentally observed increase 
in the cross section beginning at ^ 690 eV. The experimentally measured EA from n = 2 
excitations could also be systematically enhanced by field ionization, which inhibits radiative 
stabilization for excitation to n > 44, but this effect should be very small since such states 
ar e exp e cted t o autoionize even without external fields. The distorted wave calculations 



of 



Derd (120071 ) show EA of similar magnitude to what was measured ex perimentally, but 



starting at somewhat higher energies of ~ 750 eV. 



Pindzola et al. 



(119861 ) predict additional 



Tom 2 p — )■ 4 p, 2p — > Ad, and 2s — )■ 3d excitations to be found in the range 880 - 900 eV 



EA 

and iDerd ( l2007l ) also predicts an increase in the cross section due to EA at those energies. 
However, we do not see a significant increase in the experimental data in that energy range. 
A possible explanation for this is that the theory has overestimated the branching ratio of 
the intermediate state for autoionization versus radiative stabilization. 

In the energy range of 650 - 950 eV resonances are observed in the cross section, as 
shown in more detail in Fi gure El These are due to the REDA and READ! processes which 



result in single ionization (ILaGattuta &: Hahn 



1981 



Henry fc Msezane 



1982 



Miiller et al. 
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1988bl Jah. 



At collision energies near 1100 eV direct ionization of the n = 2 electrons becomes 
possible, but no corresponding feature is observed in the single ionization cross section 
at this point. This is not surprising since the excited states created by these L-shell 
ionizat ions are expected to au toionize producing a net double ionization, as is discussed 



below ( iKaastra fc Mewe 



19931). 



The experimental cross section was used to determine the plasma ionization rate 
coefficient ai(Tc). This calculation involves multiplying the ionization cross section by the 



electron-i on relative v e 
Following iFogle et al. 



ocity and integrating over the Maxwellian velocity distribution. 



( I2OO8I ) we truncate the integral over energy at Eq + Qk^Te^ where 
Eq = 330.79 eV is the ionization potential of Fe^^"'", Tg is the temperature at which ai(Te) is 
to be calculated, and fee is the Boltzmann constant. The cross section was measured only 
np to E = 2300 eV. Thus, the rate coefficient for Te < 3.8 x 10^ K could be calculated 
without extrapolation. In order to calculate ai for > 3.8 x 10^ K we assum ed that the 



Derel (120071 ) 



falloff of 0"! at high energies is the same as that given by the FAC calculation of 
and scale that cross section to our experimental results to extrapolate the measured cross 
section beyond 2300 eV. The uncertainty on the experimentally-derived rate coefficient is 
±12% due primarily to the systematic uncertainty in the ion current calibration. 

Figure [3] compares the plasma rate coe fficien t frorn the present measurements to 
those obtained from the FAC calcul ation of iDerd f|2007l ) and from the fit to the crossed 



beams experimental cross secti on (IGregorv et al. 



Arnaud fc Raymondl (119921 ) and 



Mazzotta et al. 



1998 



19871) used in the CIE calculations of 



). The ionization equilibrium depends 



most strongly on the rate coefficients near the peak ion abundance. Recent CSD calculations 
show that Fe^^"^ has a greater than 1% abundance in the temperature range from 8 . 6 x 10 ^ K 



to 2.5 X 10^ K, with a maximum ion abundance near 1.6 x 10^ K (IBryans et al 



20091). 
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The rate coefficient from the lArnaud &: Raymond (Il992l ) fit is about 30% larger than our 
results throughout this temperature range. The rate coefficient from the distorted wave 
calculation is about 15% smaller than the present measurement at the temperature of 
maximum abundance (Figure [3]). The difference is likely due to the exclusion of EA from 
the n = 3 level in the calculation, as well as potential systematic effects of field ionization 
in the experiment. Since these processes could not be distinguished experimentally, the up 
to ~ 25% difference represents an upper bound on the possible effect of EA from n = 3 
levels. It appears that including EA from the n = 3 level in the calculation of ionization 
cross sections and rate coefficients is necessary for more accurate EISI rate coefficients. 



In ord e r to p roduce a ffi to the rate coefficient we used the Burgess- TuUy type scaling 



from 



Dere 



(120071 ) and fit the scaled rate coefficient with a fifth order polynomial. The 
temperature was scaled as 



X = 1 — 



In 2 



ln(t + 2)' 

where t = k-oTe/EQ. The rate coefficient was scaled as 



(4) 



P 



(5) 



where Ei{l/t) is the ffist exponential integral. The scaled temperature x and rate coefficient 
p can be inverted to reproduce Tg and ai using 

In 2 



Eo 



exp 



and 



1 — X 



ai{%)=t-'/'E^'/'E,{l/t)p. 



(6) 



(7) 



The coefficients for a fifth order polynomial fit to p{x) are given in Table El The fit 
reproduces the experimental ai(Te) to 1% or better accuracy over the temperature range 
Tp = 2 X 10^ - 1 X 10^ K. 
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Results for Double Ionization 



The double ionization cross section in the energy range 500 - 2500 eV is shown in 
Figure |H The points show the experimental values. The dotted lines indicate the la 
systematic un certainty from the BPM calibration. The threshold for direct double ionization 



is 698.83 eV flRalchenko et al. 



20081 ). but the measured cross section is not significantly 



different from zero below the threshold for single ionization of n = 2 electrons near 1000 eV. 

The solid line in Figure H] shows the cross section for single ionization of an n = 2 
electron calculated in the distorted wave approximatio n using the GIPPER code of the 



LANL Atomic Physics Code suite ( iMagee et al. 



19951 ) scaled to reflect the probability 



of autoionization relative to radiative stabilization. The scaling factor s used were 0.99 



for io nization of a 2p electron and 0.92 for ionization of a 2s electron (IKaastra fc Mewe 



19931 ). Thus, the solid line represents the cross section for double ionization proceeding 
through EISI of an L-shell electron followed by autoionization when the "hole" is fllled 
in. The ElDl cross section is clearly dominated by L-shell single ionization followed by 
autoionization, as seen by the good agreement between the experimental cross section and 
the calculated cross section for this process. Such indirect processes are known to become 
i ncreasingly importan t relative to direct i onization as the ch arge state of the ion increases 



(IMiiller fc Frodl 



1980; 



Miiller et al. 



1985 



Stenke et al. 



19991 ). 



7. Summary 



We have measured the EISI and EIDl cross sections of Fe^^^. In the single ionization 
case we found that the cros s section was ^ 30% smaller than that measured in a single-pass. 



crossed beams experiment ( IGregory et al. 



19871 1. The difference is believed to be due to the 



presence of metastable ions in the single-pass experiment and the absence in the current 
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experimen t. The sing le ionization cross section is in agreement with the theoretical cross 



section of 



Derd (120071 ) to within about ±20%. The remaining differences are probably 



due to differences in EA between the calculation and the experiment. At low energies the 
calculation leaves out EA from excitation of 3s electrons. At higher energies the calculation 
seems to overestimate the EA contribution. Our results show that the Maxwellian EISI rate 
coefficient from calculations of 



Derd (120071 ) may be underestimated by as much as ~ 25% 
over the temperatures at which Fe^^"*" is abundant, partially due to the need to include EA 
from the n = 3 level in the calculations. The experimental results for double ionization 
show that the cross section is dominated by single ionization of an L-shell electron which is 
followed by autoionization. This is consistent with previous results for lower charge state 
ions, which showed these processes tend to become more important as the ion charge state 
increases. 
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A. Correction for EII at the Reference Energy 

Equation (jj]) includes a correction term that accounts for a non-zero EII cross section 
at the reference energy. This term includes a relativistic factor that has been omitted in 
previous work. This appendix briefly explains the origin of this term. 

The ionization count rate at the measurement energy is given by 

RT{E^) = (ait;,ei) (i?nO<Ar-L[l - (3-MEm)] + Res- (Al) 
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Here Res represents the background count rate due to electron stripping off the rest gas. 
The other symbols are defined in the discussion following equation ([2]). For simplicity we 
ignore pressure fluctuations and the decay of the ion beam and assume that i?Es and A^i 
are the same at the measurement and reference energies. Systematic errors ar ising from 



this a ssumption can be corrected for separately using the method described in 



Hahn et al. 



(120101 ). With these simplifications the background count rate measured at the reference 



energy is given by 

R'^iE^) = (orit;,ei) {E,)nlNrL[l - A/3c(^r)] + Res- (A2) 
Subtracting equation (IA2p from equation flAll) and solving for (aifrci) (-Em) gives 

We recover equation ([2]) after dividing by the relative velocity. In terms of v/c, the velocities 
in the present experiment were small and [1 — PiPe{E)] ^ 0.99. Hence this additional 
relativistic correction has a nearly negligible effect. 



B. List of Abbreviations 

BPM Beam Profile Monitor 

CEM Channel Electron Multiplier 

CIE Collisional Ionization Equilibrium 

CSD Charge State Distribution 

DES Double Electron Stripping 

DR Dielectronic Recombination 

DW Distorted Wave 

EA Excitation Autoionization 

EIDI Electron Impact Double Ionization 
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EII Electron Impact Ionization 

EISI Electron Impact Single Ionization 

REDA Resonanat Excitation Double Autoionization 

READI Resonant Excitation Auto-Double-Ionization 

SEC Single Electron Capture 

SES Single Electron Stripping 
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Table 1:: Sources of Uncertainty. 





ll/oLlilidLcU. XC liL-Ci LdliiL V 


Ion current measurement 


12% 


Electron density 


1% 


Detection efficiency^ 


< 0.1% 


Interaction length spread^ 


< 0.1% 


Counting statistics 


1% - 10% 


Background correction^ 


0-3% 


Quadrature sum 


12% - 16% 



^Although negligible, these un certainties are listed here for comparison with the similar 



experiment of 



Hahnetal 



(|2010[ ) where the expansion of the ion beam introduced small 
additional uncertainties. In the present experiment continuous cooling of the ion beam 
inhibited beam expansion making these measurements more precise. 

^The 3% uncertainty applies to data for E > 950 eV for single ionization and E > 1400 eV 
for double ionization where the background correction could not be applied. 
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Table 2:: Fifth-order Polynomial Fitting Parameters to 
Reproduce the Scaled Ionization Rate Coefficient p(x) 
(see equations H] and [5]) . 



i ai 

19.4438 

1 -76.2632 

2 376.369 

3 -819.161 

4 827.141 

5 -322.187 

1=0 
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T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




j—i± I I I I I I I I I I I I I I I I I I L 

500 1000 1500 2000 

Energy (eV) 



Fig. 1. — The EISI cross section for Fe^^+ forming Fe^^"^ is shown here. The filled circles 
indicate the experimental values and the error bars at selected points illustrate the Icr sta- 
tistical uncertainty. The dotted curves show the la range of systematic uncertainty from 
the calib ration of the ion cur rent measurement. The diamonds display the experimental 
results of 



of 
of 



Gregorv et ah 



(119871) and a fit to that data that was used in the CIE calculations 



Arnaud fc Ravmondl (119921 ) is indicated by a dashed line. The distorted wave calculation 



Derd (120071 ) using the FAC code is denoted by the solid curve. 
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5 
650 



This experiment 
FAC (Dere 2007) 



700 750 800 850 
Energy (eV) 



900 



950 



Fig. 2. — This figure shows the EISI cross section of Fe^^^ focussing on the resonances in 
the range 650 - 950 eV. The error bars indicate the la statistical errors, which are small 
co mpar e d to t he size of the resonances in this energy range. The distorted wave calculation 



of 



Derd (120071 ) is also shown. 
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Fig. 3. — The EISI rate coefficient for Fe^^+ forming Fe^^+ derived from the experimental 
data (sohd hne). A lso shown are rate coefficients derived from the earlier experiment of 



Gregory et al. 



(119871 ) by Arnaud & Raymond (1992; dashed line) and from the FAC calcula- 
tion by Dere (2007; dash-dotted line). The relative differences between these and the present 
experimental rate coefficients are also illustrated. Arrows are placed next to the curves to 
clarify the axis from which the values are to be read off. The vertical dotted lines show the 
temperature range where the equilibrium ion abundance of Fe^^"*" is > 1% with the m iddle 



dotted line showing the temperature of maximum ion abundance (IBryans et al 



20091 ). The 



systematic error on the experimentally derived rate coefficient is about ±12%. 
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1.5 




500 1000 1500 2000 2500 

Energy (eV) 



Fig. 4. — The EIDI cross section for Fe^^"'' forming Fe^^^ is shown here. The la statistical 
uncertainties are indicated by the error bars. The dotted curves illustrate the systematic 
la uncertainty from the ion curr ent measurement. Also shown is a calculation using the 



GIPPER code flMaeee et al. 



19951 ) for the double ionization cross section due solely to single 



ionization from an L-shell electron followed by autoionization. 
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